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Confronting pNGB quintessence with data 
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We analyze the observational constraints on the model where a pseudo-Nambu-Goldstone boson 
(pNGB) plays the role of dark energy. The constraints are derived by using the latest Gold set of 182 
type la supernovae and the CMB shift parameter. We allow for both the initial value of the scalar 
field and the present value of the energy density in the pNGB to vary. We find that - compared to 
previous analyses - the allowed portion of parameter space has shrunk around the region where the 
pNGB does not evolve significantly. 
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I. INTRODUCTION 

Since the type la supernova observations of [l], 0, Q , 
an intense activity has been devoted to the search of an 
explanation of the accelerated expansion of the Universe. 
If gravity is described by Einstein's General Relativity 
and the effects of inhomogeneities can be neglected, then 
acceleration must be due to a dark energy component 
that represents roughly 70% of the matter content of the 
Universe. Current data tell that the equation of state 
parameter w = p/p of dark energy has to obey w < — 
0.7 0. 

The simplest explanation of cosmic acceleration is a 
cosmological constant, the energy of vacuum, with mag- 
nitude 



A ~ (2 x 10~ 3 eV) 4 



(1) 



and equation of state parameter w = — 1. This solution 
is attractive in many respects, both for its simplicity (a 
single parameter is needed to describe it) and for its ex- 
cellent agreement with data. It is however hard to justify 
from a theoretical standpoint. Quantum fluctuations of 
matter, indeed, give contributions to the vacuum energy, 
and very precise cancellations are needed to keep this en- 
ergy at small values. Unfortunately the Standard Model 
does not display any of these cancellations at least up to 
the scales that have been probed in collider experiments, 
about 60 orders of magnitude beyond the value of ([T]). 

For this reason, soon after the release of [l|, 0, [ij] , people 
have started to look for alternative scenarios, and a wide 
interest in quintessence models has emerged [5[ (for a re- 
cent review, see Q). The philosophy behind quintessence 
is the following. First, it is assumed that some mecha- 
nism is able to fix the energy of the ground state of the 



Universe to zero 1 . Then the existence of a new degree 
of freedom (quintessence) is postulated: quintessence is 
supposed not to have yet relaxed to its vacuum, so that 
its energy density is responsible for cosmic acceleration. 
Quintessence has w ^ —1 as a distinctive prediction, and 
is usually described by some scalar degree of freedom <f) 
endowed with some potential V ((f)). V ((f)) has to be very 
flat, if we want w to be sufficiently negative. 

It is possible to write down a virtually infinite number 
of quintessence potentials V (4>). However, only for few 
of them the flatness of the potential is not spoiled by ra- 
diative corrections and the exchange of quanta of (f> do 
not give rise to an (unobserved) fifth force Q • Those few 
potentials are more motivated from a theoretical point 
of view than the others. This is especially true for the 
pseudo-Nambu-Goldstone boson (pNGB) potential of Q 
that has all the good qualities of radiative stability that 
anybody who believes in quantum mechanics might re- 
quire. 

In the present paper we analyze the parameter space 
of pNGB quintessence Q in the light of the most re- 
cent observations, in particular those from supernovae. 
Our approach is orthogonal to the "model independent" 
approach recently taken on the subject by many investi- 
gators (see for instance Q), and is admittedly based on 
a theoretical prejudice in favor of radiatively stable po- 
tentials. To our knowledge, the most recent complete 
anal yses of thi s model date back to about five years 
I. Given the recent developments of the 
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observational situation, we believe that it is important 
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1 It is often stated that finding a mechanism that fixes to cosmolog- 
ical constant to zero should be easier than finding a mechanism 
that fixes it to some very small nonvanishing value. Let us note 
here that this is not what usually occurs in Quantum Field The- 
ory: if it is possible to find a symmetry that fixes some quantity 
to zero, it is typically straightforward to break such a symmetry 
so that this quantity can be kept small in a controlled way. 
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to perform an analysis of the model in which the latest 
data are taken into account. 

In the next section we briefly describe the properties 
of the model of pNGB quintessence. Then in section III 
we present the observational constraints from supernovae 
and from the CMB shift parameter. In section IV we 
discuss our results before concluding in section V. 



II. THE PNGB POTENTIAL 

The use of pNBGs has been first proposed in order to 
realize a technically natural model of inflation in 13[ and 
has been subsequently considered for dark energy in 
The model is characterized by a pseudoscalar field <j) with 
a potential that can be well approximated by 



V(cf>)=(x' 1 [cob fa//) 



1] 



(2) 



where we have neglected the contributions by higher 
harmonics (this is supposed to be a good approxima- 
tion as long as / is sufficiently smaller than the Planck 
mass [lj]). The potential is generated by the breaking of 
a shift symmetry <f> — > </>+constant, and for this reason it 
is radiatively stable. 

The cosmological evolution of this model is in general 
determined by four parameters: the quantities /i and / 
and the initial conditions 4>i n and <j)\ n . Due to the high 
expansion rate of the Early Universe, we assume <j>- in = 0. 
One more free parameter is eliminated if we insist that 
fi^ (the current ratio of the amount of dark energy over 
critical energy) is equal to 0.7. As a consequence, if we 
assume = 0.7, the model is described by only two 
parameters that can be taken to be / and <j>- iri . A detailed 
analysis of the dynamics of the pNGB zero mode can be 
found in 



111 ]. Due to periodicity of the potential 



takes values between and 2irf. In addition, if we take 
into account the indication from String Theory 14( that 
/ cannot be larger than Mp ~ 2.4 x 10 18 GeV, then the 
parameter space of the potential turns out to be compact. 
This implies that, at least in principle, all of it can be 
excluded by observation, and that pNGB quintessence 
in its simplest version can be ruled out. For this reason 
we find this model even more attractive (although there 
are ways to evade the constraint / < Mp [15|), and we 
believe this is an additional motivation for studying it in 
detail. Also, for this reason we will restrict our study to 
the region / < Mp. 



III. OBSERVATIONAL CONSTRAINTS 

We consider the quintessence field <j> with potential (J2J) 
in a flat Friedmann-Robertson- Walker Universe. The 
equations of motion are given by 



(3) 



where H — a/ a is the Hubble constant, p m is the energy 
density of nonrelativistic matter and, since we will be 
dealing only with the dynamics of the late Universe, we 
neglect contributions by radiation. We solve these equa- 
tions numerically to find the evolution of the scale factor 
function of time. 



A. Type la supernovae 

First, we investigate constraints from the observation 
of type la supernova from the data set [l6j |. which is a 
compilation of old data [3| byHZS team, first year Su- 
perNova Legacy Survey data [l7j and recent observations 
of 21 new supernovae [4}. For our analysis we will con- 
sider only the 182 "high confidence" Gold SN data with 
z > 0.0233. Although most of the SNe have z < 1, there 
are 16 SNe with z > 1. This is the most up to date su- 
pernova data set available in the literature. This data set 
has been recently used in 3, 1^| to study the observa- 
tional constraints on different parametrizations of dark 
energy. 

For a particular cosmological model with parameters s 
the predicted distance moduli are given by 



Ho(z, s) = m - M = 5 log 10 (j^f) + 25, (4) 

where m and M are the apparent and absolute magnitude 
of distant supernovae. d,L is the luminosity distance given 
by 



d L {z) = {l + z) 



dz 



H(z') 



(5) 



and depends only on the expansion history of the Uni- 
verse from redshift z to today. Assuming that all the dis- 
tance moduli are independent and normally distributed 
the likelihood function can be calculated from the chi- 
square statistics C cx exp(— x 2 /2), where 



, 4>in, /, Ml H ) 



182 l,,obs 



Moi) 2 



(6) 



'Qi 
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Here /zj^ s and aoi are the measured value of the distance 
modulus and the corresponding uncertainty for the i— th 
supernova. /iQi S and am, as well as the redshift Zj are 



found from the data set fiff is calculated by using 
eq. ([5]), where H(z) is obtained by numerically solving 
the background evolution equation ([3]). We marginalize 
the likelihood over the nuisance parameter H$ |2j]. 

As we have noted earlier, the model has four param- 
eters: /, /i, 4>i n and </>;„. We assume 4n n — and allow 
the system to evolve until Sl^ = 0.7 today. This leaves us 
with two parameters and we choose them to be / and 0; n . 
We plot the resulting confidence contour in figure [1] The 
upper left portion of the plot corresponds to the part of 
parameter space where fi^ does not reach the value 0.7. 
In this part of the parameter space, the scalar field rolls 
quickly to the minimum and oscillates around it, behav- 
ing like matter. 




FIG. 1: The shaded areas at the bottom right of the figure 
denote the 1 a and 2 a confidence level regions for £1^ = 0.7. 
The upper left part of the plot corresponds to parameters 
for which fi^, never reaches 0.7. The white dashed line corre- 
sponds to the value of the parameters for which wq = —0.965. 

The dark areas at the bottom right part of the plot give 
the 68.3% (la-) and the 95.4% (2<r) confidence level re- 
gions. The 3 a contour runs between the 2 a contour and 
the boundary of the forbidden region (we do not show it 
in the plot for clarity). As we go closer to the boundary 
of the forbidden region, the value of \ 2 increases sharply. 
Even if the 95.4% confidence level area seems to cover 
almost all of the allowed region, this is actually not the 
case. Indeed, there is a part of the parameter space, be- 
low the boundary of the forbidden region, where Q^, goes 
across 0.7 several times, as an effect of the oscillations of 
cj). This means that points with the same values of <f>- m , f 



and ilef, can correspond to different histories, depending 
on the number of times has gone across the value 0.7. 
We have computed the value of \ 2 m the case where 
has crossed VL^ = 0.7 more than once and we have found 
that this part of the parameter space is excluded at more 
than the 3 a level. 

Once fi,£ = 0.7 is fixed, the constraints on the param- 
eters 0i n and / can be converted into constraints on the 
plane (/, \x). We plot the 1, 2, and 3u contours in fig- 
ure [3 
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FIG. 2: 1, 2, and 3 a constraints on the plane (/, n) for fi^, = 
0.7. The lower part of the plot corresponds to values of the 
parameters for which fi^ cannot reach the value 0.7. 

The thick white dashed line in the bottom right part 
of figure [T] corresponds the part of parameter space that 
gives Wq = —0.965 (where wq is the current value of 
the equation of state parameter). According to [2o| . 
w < —0.965 is the most optimistic constraint (at the 
95.4% level) that we might obtain from future observa- 
tions, should they converge to the regime where dark en- 
ergy shows no evolution. Therefore the dashed line in fig- 
ure [T] gives the most stringent 2 a constraint that we can 
expect to put on the parameters of pNGB quintessence. 

We have also considered the case where the value of 
0^ is allowed to vary. In this case we have fixed f = Mp 
while keeping </>; n variable. In the left panel of figure [3] 
we show the 1, 2, and 3er contours related to supernova 
observations on the (O^, <j>- m ) plane. The shaded upper 
right part of the plot is excluded since the corresponding 
value of n<£ cannot be reached. The contours are essen- 
tially vertical and centered around ~ 0.67. However, 
at larger values of f2^, somehow larger values of 4>- m (cor- 
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responding to some evolution in the quintessence field) 
are allowed. The best fit is at = 0.67, </>;„ = (so 
that the pNGB sits at the top of its potential and be- 
haves as a cosmological constant) with x 2 = 159.6 for 
180 degrees of freedom 2 . For smaller values of / (figure 
not shown), the contours have the same shape, although 
they shrink along the 4>[ n direction. 
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FIG. 3: Constraints on the pNGB parameter space for / = 
Mp. Left panel: 1, 2 and 3<j confidence level contours from 
supernovae only (dashed lines) and from the CMB shift pa- 
rameter © only (solid lines). Right panel: 1, 2 and 3cr con- 
tours from the joint analysis. 



the left panel of figure [3] we plot the resulting confidence 
contour arising from this constraint. Note that the CMB 
contours favor a value of that is slightly larger than 
that favored by the SN data. 

Once the SN and the CMB constraints are combined, 
we obtain the plot shown on the right panel of figure [3] 
Since we are assuming a flat Universe, imposing the shift 
parameter constraint does not reduce significantly the 
area of the allowed region, and indeed the SN and CMB 
constraints are not orthogonal. However, the CMB con- 
straint helps eliminate the part of parameter space at 
large and large 4>- ln that is available at the 2er and 
3cr level if only SN constraints are taken into account. 
Moreover, once the CMB constraint is added, the best 
fit point is not any more at = 0, but at the point 
(j)- m = 1.25 Mp,ft = 0.71 where x 2 = 161.6. This should 
be compared to \ 2 — 162.9 found at Q,<p — 0.71 when the 
constraint <f>- in = is imposed. Since 4>m implies 
that (j> is rolling, the combination of CMB and SN data 
seems to hint at some evolution in dark energy. How- 
ever, this hint should be taken with a grain of salt, since 
it emerges when we join two data sets that are not ex- 
actly compatible, as shown by the increase of ~ 2 units 
in x 2 when we add the single CMB point to the SN data 
set (as stated above, the best fit point for SN data only 
has x 2 = 159.6). 

IV. DISCUSSION 



B. CMB shift parameter 

In addition to the SN data, we use Cosmic Microwave 
Background (CMB) data to constrain the model. In par- 
ticular we derive constraints from the CMB shift param- 
eter 1Z, that measures the shift in the angular size of the 
acoustic peaks of CMB when parameters of the theory 
are varied. 1Z is independent on the present value of the 
Hubble constant, and is given by 



K 



(7) 



where z cm b is the redshift of recombination. By using 
WMAP 3rd year data the value of the shift parameter has 
been extracted as 1Z — 1.70±0.03 for z cm b = 1089 [2l|. In 



2 The best-fit x 2 for the older dataset [J was of \ 2 = 178.1 for 155 
degrees of freedom. The lower value of \ 2 for the current data 
set can be largely attributed to more conservative assumptions 
on the dispersions coi- 



The data presented above indicate that, if we use su- 
pernova constraints only, the parameters that yield the 
best fit to data are those where the field cf> sits at the 
top of the cosine potential, thus mimicking a cosmologi- 
cal constant. If we use also the constraint from the CMB 
shift parameter, a slowly rolling pNGB is slightly pre- 
ferred to a constant one. 

In order to see how the new data of Q improve the 
constraints on the model, we can compare our results 
with those of previous analyses. In Waga and Frie- 
man have studied the constraints from the 1998 super- 
nova data of Riess et al. [1| , together with the statistics of 
gravitationally lensed quasars. Comparison of our results 



with those of lOj is complicated by different assumptions 
on the parametrization of the model. Indeed, in [lCj the 
value of 4>in is fixed to 1.5 /, so that £1$ is not a free pa- 
rameter, but is function of / and /i. Nevertheless, some 
comparison is possible: in the parameter space of 10 1 
there is still room at 2 a for a small region where the 
scale factor of the Universe is currently decelerating. In 
our analysis (see figure 2]) this is not possible any more 
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at the 2 a level, even if it is still allowed at 3 a. 

In a wider portion of the parameter space is an- 
alyzed (and a different data set [2| is used), that shows 
that a part of parameter space where <p has performed 
half oscillation is allowed at the 2 a level. As we have 
stressed in the section III, the current data do not allow 
for this possibility any more. 



In 



12j, a detailed study of the parameter space 



of the model has been performed by taking into 
account the constraints from CMB observations of 
BOOMERanG [22] and MAXIMA Q. Constraints on 



this quintessence model were derived from its effects on 
integrated Sachs- Wolfe effect as well as from its effects 
on the location of the first peak. In figure [5] we show 
the 1, 2, and 3 a constraints in the (/, /i) plane obtained 
by our analysis. Comparison with figure 5 of shows 
that the more recent data improve by a factor of three 
or so the constraints on \x. 

Let us also note that a quintessence field that is climb- 
ing up the potential could mimic w < — 1 [24|, 25 1 and pos- 
sibly offer a better fit to data (see however 26|). In our 
case this is possible only if has already gone through 
a maximum. As we have mentioned above, this case ap- 
pears to be not realized at the 3 a confidence level for 
a cosine potential. Indeed, a more asymmetric potential 
(such as that pictured in figure 1 of [25j) is needed to im- 
part a sufficiently large velocity to the field and improve 
the fit to data. 
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FIG. 4: 1, 2, and 3 a contours in the plane (fio.E, wo) for the 
pNGB model with / = Alp (dashed, thicker lines) and for 
dark energy with constant equation of state (solid lines). 

Finally, it is instructive to discuss the current value of 
the equation of state parameter wq as obtained in the 
pNGB model and to compare it with the value of Wq ob- 
tained by assuming that it is constant throughout the 



evolution of the Universe. In figure [4] we show the plots 
of the 1, 2 and 3 a confidence level contours in the plane 
(fi^, wo) both for the case of pNGB quintessence with 
/ = Mp and for a model with constant wo (only super- 
nova data are used to compute the contours in figured]). 
In the case of constant wq the contours are quite tilted, 
allowing for a value of wq significantly different from —I 
only if Que gets very close to unity. Indeed, if dark en- 
ergy has a constant equation of state parameter different 
from —1, then a larger amount of dark energy is needed 
to get the same averaged value of w. In the case of a 
pNGB, however, the value of Wq can be significantly dif- 
ferent from —1 even if w was close to —1 in the past. 
As a consequence, one can obtain the required averaged 
value of w even without requiring that is extremely 
close to unity. 

V. CONCLUSIONS 



We have analyzed the portion of parameter space avail- 
able for the model of pNGB quintessence. Our work ex- 
tends the previous studies on the subject by allowing 
both for variations in the initial value of the zero mode 
of the pNGB and for variations in the current value of 
£1^. Using the most up-to-date supernova data, we have 
shown that the parameter space of the pNBG poten- 
tial is significantly constrained around the region where 
quintessence is sitting on the top of the cosine potential 
or slowly evolving along it. At the 95.4% level, previ- 
ous analyses on the subject [13, U| were still allowing 
the current value of w to be larger than —1/3 or even 
the possibility that quintessence had already performed 
a half oscillation about its minimum. Current data do 
not allow this any more. 

We have also observed that, when CMB and SN con- 
straints are joined, an evolving pNGB provides a slightly 
better fit to data than a pNGB stuck at the top of its 
potential. 

Let us finally discuss future perspectives. Already now, 
data tell that / cannot be smaller than about a third of 
Planck mass (unless we fine tune <pi n to be very close to 
the top of the potential). As shown in figure [TJ future 
data might constrain / >, Mp/2, leading to some tension 
with the requirement / ^$ Mp from String Theory 14 1. 
One might wonder if this will be enough to consider the 
model "finely tuned", and to start to consider alterna- 
tive options Q as more natural. But one can also take 
a more optimistic approach: maybe future data will show 
that cosmic acceleration is sourced by an evolving, radia- 
tively stable pseudo-Nambu-Goldstone boson. 
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